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Abstract 

We present an analytic study of the phase diagram of the one-dimensional 
t — J model and a couple of its cousins. To deal with the interactions in- 
duced by the no double occupancy constraints, we introduce a deformation 
of the Hubbard operators. When the deformation parameter A is small, 
the induced interactions are softened, accessible by perturbation theory. We 
combine bososnization with renormalization group techniques to map out the 
phase diagram of the system. We argue that when A ^ 1 , there is no essen- 
tial change in the phase diagram. Comparison with the existing results in the 
literature obtained by other methods justifies our deformation approach. 
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I. INTRODUCTION 



A. Historical review 

Since the discoveries of quantum Hall effects and high Tc oxides in 1980's, the strongly 
correlated systems have been of great interests both theoretically and experimentally. As 
far as the high problem is concerned, the t — J model is believed to be the appropriate 
starting model Hamiltonian, because it captures the essence of the interplay between charge 
and spin degrees of freedom in superconducting Cu oxidesS'i. Although high-Tc cuprates 
are (at least) two-dimensional systems, it is very interesting to study its one-dimensional 
(ID) counterpart. As argued by Anderson^, two-dimensional strongly correlated systems 
may share some properties of ID systems. In addition, the physical understanding of the 
ID systems is also extremely helpful for the study of the ladder systems, which have been 
realized experimentally and have attracted a lot of attention in recent yearsiHl. 

In the ID systems, the phase space of particle scattering is highly restricted. The oc- 
currence of a single scattering event will spread quickly among all other particles, which 
invalidates the concept of individual excitations. Consequently, we are often confronted 
with correlated collective excitations. On the other hand, in some cases we can benefit from 
such phase space restriction. Namely, the many-particle scattering matrix could be nicely 
decomposed into the product of two-particle ones which satisfy the so-called Yang-Baxter 
integrable conditionsi. This property provides us with the possibility to exactly solve some 
ID models, e.g., the Hubbard mo del0, the Heisenberg model&i, and the supersymmetric 
t — J model0. The exact solutions in turn provide us with powerful guidelines to develop 
and to justify certain approximate schemes for other problems. 

In some sense, the ID t — J model could be viewed as a descendant of the Hubbard model 
in large on-site repulsion limit. Namely, the strong coupling limit of the Hubbard model 
can be mapped into the weak coupling limit of the t — J model. Naively, one may speculate 
that the integrability of the ID Hubbard model would be inherited by the t — J model in 
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the whole parameter space. Unfortunately, this speculation is not correct: The t — J model 
is only integrable at two special points in parameter space. The reason for this difference 
is that, in contrast to other ID integrable models, the Hilbert space of the t — J model is 
highly constrained: Double occupancy of any site is completely excluded. Furthermore, the 
integrable points of the t — J model are located in the strong coupling regime (| J| = 2t in 
our convention below), not in the weak coupling limit. Therefore, the integrability of the 
supersymmetric t — J model is not simply inherited from the Hubbard model. Rather it is 
better to be viewed as a separate miracle of the interacting ID many-particle system. Since 
the ID t — J model can not be exactly solved at a generic point in parameter space, the 
analytical studies of the t — J model have been a painstaking task even in the ID case. 

To illustrate the points more clearly, let us take a close look at the t — J model. The 
model delineates the behavior of hard core fermions on a discrete lattice, and the dynamics 
is given by the model Hamiltonian 

Htj = + h.c.)V + jE ■ S,+i. (1) 

Here V is the projection operator that prohibits double occupancy of any site, a and a = —a 
the spin orientations (with cr=l for |, and -1 for |); t is the hopping amplitude and J the 
anti-ferromagnetic (J > 0) or ferromagnetic (J < 0) coupling. Due to the aforementioned 
constraints, at each site the states \a > can only be one of the following three possible states: 
with a =T; i, and a = (empty). This Hilbert space is neither fermionic nor bosonic. One 
can check that the projection operators x"^ = |a >< 6| close, under commutation and anti- 
commutation, to form a semi-simple supersymmetric Lie algebra, the Spl{l, 2) given by the 
relations^ 

{xf , Xf}± = S.^ixfS'^ ± 5'^'), (2) 

where Xi^ x^'^ ai'G fermionic operators that, respectively, create and annihilate a single 
electron. The bosonic operator are identified as the generators of the group SU{2). 
Using these operators, the t — J model can be neatly written as 
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HtJ = -t ixTXjl + h.c.) + J E xf'x^i, (3) 

i,o-=t,J, j,(T,a' 

in terms of the bihnears in the generators of Spl{l,2). But the price we have to pay is 
to introduce both fermionic and bosonic operators simultaneously. As a consequence, it is 
difficult to make a simple, controlled approximation in this representation. To overcome the 
difficulties associated with the no-double-occupancy constraints, the slave boson and slave 
fermion methodili and, more recently, the supersymmetric Hubbard operator methodlil 
have been invented to treat the t — J model, with the hope of the mean field ground state 
being relevant to the high-Tc problem. However, after more than one decade effort, it seems 
that a reliable ground state is still elusive. 

B. Deformed Hubbard operators 
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The seminal work by Jordan and Wigne and, later, by Lieb, Schultz, and Matti 
provides an alternative idea to handle the above hybridized situation in statistics: Namely 
the spin operators are uniformly expressed in terms of fermionic operators, though the spin 
systems are neither bosonic nor fermionic ones. In the same spirit, one can also rewrite the 
t — J model in terms of fermions exclusively. 

In addition to rewriting the magnetic interactions using the fermionc realization of the 
spin operators 

= \c]^^o.pCjp, (4) 
one may introduce as well the Hubbard operators^ 

c]a = 4^(1 - nja) (5) 

and rewrite the hopping terms in terms of them. These operators also realize the constraints 
that exclude double occupancy on each lattice site. In this way, one gets a formulation of 
the t — J model completely in terms of fermionic operators. 
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However, one immediately sees that the old hopping terms will induce extra four-fermion 
and six-fermion interactions. These interactions are "hard" ones, in the sense that their 
strengths are exactly the same as the hopping amplitude t. This fact defies the attempts to 
treat the additional four- and six-fermion terms perturbatively. Therefore, at first glance, 
it looks silly to adopt this strategy to solve the t — J model unless new techniques can be 
invented to make the induced interactions tractable. 

In the present paper we propose a technique that allows us to deal with these induced 
four- and six-fermion interactions. The key point is to use the idea of " adiabatic continuit'ij'' 
to soften the above-mentioned interactions induced by the no-double-occupancy constraints. 
Namely, we propose to introduce the following deformed Hubbard operators 

cl = clil - An,,) (6) 

with a deformation parameter < A < 1. When A approaches to unity, we recover the 
genuine Hubbard operators (^. For < A < 1, there is a non-zero probability to allow 
leakage into states with double occupancy. With these deformed Hubbard operators 
replacing the genuine Hubbard operators (|^) in the hopping terms, we get a deformation 
of the original t — J model. The deformed model has the advantage that for small A, 
the induced four- and six-fermion interactions are no longer "hard". This is because these 
interactions have strengths proportional to the deformation parameter A and, therefore, is 
tractable in the sense of perturbation theory when A is small. 

Though small values of A may not be "physical" , after extracting possible structures in 
the phase diagram for small A, we analytically continue our results back to A = 1. The 
fundamental assumption underlying this continuation is the adiabatic continuity, namely 
that when the Hamiltonian of the model is adiabatically changed with A varying from a 
small positive value to unity, there is no essential, qualitative change in the phase diagram of 
the model, though various phase boundaries in parameter space may undergo a continuous 
deformation. Historically, our idea of considering a deformed model is parallel to the ideas 
that underlie the replica method in treating disordered system, or the large expansion in 
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field theory. Actually, even in the field of ID exactly solvable models one can find a precedent: 
Yang and Yang|l proposed the XX Z model as a deformation of the XXX model, i.e. the 
spin-i ID Heisenberg model, and used it to justify the Bethe Ansatz method in the latter by 
first studying the large anisotropic limit and then continuing back to the isotropic limit. In 
this paper we will first discuss some simpler cases and give arguments to justify the adiabatic 
continuity assumption together with our deformed Hubbard operators. 

Of course, practically the justification may depend on how we treat the deformed model, 
which is a fully fermionized model containing four- and six-fermion interactions. In the 
present paper, we are going to combine the bosization method and perturbative renormal- 
ization group (RG) techniques to deal with the deformed t — J model. Namely we first 
bosonize the deformed model, and then use the RG fiows to map out the phase diagram of 
the bosonized model. We will argue that the phase diagram obtained in this way does not 
change in an essential way, when the deformation parameter A varies from a small positive 
value to unity. 

For convenience, we will start with a simplified model. Namely, we will first con- 
sider a model in which the magnetic spin-spin interactions are Ising-like, i.e. of the form 
JzJ^j SjSj^i- This model, together with the usual hopping term, we call the t — model. 
Then with a bit more complication, we would like to modify the isotropic magnetic interac- 
tions in eq. (|l]) to anisotropic XXZ-type interactions: 

H = -tV E(4c.+i,'^ + ^-c-)^ + J± E(^J^j+i + S]Sl,) + J, E S^S^^,. (7) 
J> j j 

This model we call as the t — J± — Jz model. The phase diagram of the SU{2) invariant 

t — J model can be obtained in the double limit with Jj_ —>■ (the isotropic limit) and with 

the deformation parameter A ^ 1 (the physical limit with no double occupancy). 

The paper is organized as follows: In section II, we discuss the phase structure of the 

extremely anisotropic limit of the t—J± — Jz model, namely, the t—Jz model. The convention 

of our bosonization scheme is also presented in detail in this section. Then the discussions of 

the phase diagram for the ID t — J± — Jz model are presented in the section III. In section IV, 
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we compare our results with other work. The discussions and conclusions are summarized 
in the section V. 



II. AN EXTREMELY ANISOTROPIC LIMIT: THE T-Jz MODEL 

A. The model 

The ID t — Jz model represents a strongly anisotropic limit of the SU{2) t — J model, in 
which only has the Ising part of the magnetic interactions been included. Without hopping, 
this simplification is significant for understanding purely magnetic interactions. However, 
with hopping the model is more interesting in that it has incorporated the interplay between 
hopping and the exchange interactions, which makes the physics of the model highly non- 
trivial. Therefore, the model has recently attracted a lot of interests^. It is known from 
the numerical studies that the low-energy physics in both the t — and t — J mo delH 
shares some common features even in two dimensions. In the real world, the possible origin 
of exchange anisotropy is the spin-orbital couplin^2l. In the extremely anisotropic limit, the 
Hamiltonian (for the t — model) reads 



'Htj^ = -tY.{clc,+,,^ + H.c.) + JzYl S]S]^, (8) 

JO" j 

= Ho{t) + U{Jz). 
Following Eq. (H), we use the representation of Sj given by 

S-=linn-nn)- (9) 

Note the appearance of the Hubbard operators in the hopping terms. It is the presence 
of the second term in eq. (§) that realizes the no double occupancy constraints. As a 
consequence, the term Ho{t) is no longer a simple hopping of fermions: More interaction 
terms with four or six fermions are induced, with strengths of the same order of magnitude 
as the hopping amplitude t. How to deal with these interaction terms is an important issue. 
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To reduce the strengths of the interaction terms induced by the no double occupancy 
constraints, we propose to deform the model Hamiltonian by replacing the Hubbard 
operators with the deformed Hubbard operators resulting in 

Ho{t) = Hh + Hi + H2 + ifa, (10) 

The Hamiltonians Hi{i = h, 1,2,3), in terms of the genuine fermion operators Cja and c]^., 
are given by 

H, = -tY.{clcj+^,^ + H.c.), (11) 
which represents the genuine hopping term, and 

Hi = tA Y,{c]aCj+i,anj+i,a + H.c.)] (12) 

-jo- 

H2 = tA J2{c]aCj+i,anj^ + H.c.)- (13) 

jo- 

H3 = -tA^ Y^{c]„Cj+i^„nj^nj+i^^ + H.c). (14) 
jo 

Here Hi and H2 are the induced four fermion repulsive interaction to prevent double 
occupancy of the same lattice site, and the H^ term is attractive, representing the effects 
from the six fermion interactions that compensate to the excessive repulsion in Hi and H2. 
It is easy to see that now in the deformed model, all the induced terms i^i, H2 and H^ 
are proportional to the deformation parameter A in Eq. (|^). If A is small, the induced 
interactions are "softened", becoming tractable in perturbation theory. In the limit of 
A — i> 1, the total effects of the three terms precisely prevent double occupancy for each 
lattice site. 

By using Eq. (|]), the exchange term U{Jz) is given by 

U{J,) = Y^(njt - nji){nj+i^^ - nj+i,i). (15) 
j 

In this way, we rewrite the t — Jz model in terms of fermion creation and annihilation 
operators exclusively. To look for the low energy effective Hamiltonian, we perform the 
standard procedure to bosonize the t — model in the next subsection. 
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B. Bosonization 



The hopping term is easily diagonahzed by Fourier transform, the energy spectrum is 
given by 

£{k) = -2tcos{ka), (16) 

where a is lattice spacing. In the ground state, all the states with momentum lower than 
the Fermi momentum kp are filled. For a generic filling factor u = N/M with the particle 
number and M the number of lattice sites, the Fermi momentum is 

TT 

kpa = -V. (17) 

To get the low energy effective action for the excitations, we only need to focus on momenta 
close to ±.kp and linearize the spectrum as 

£(±A;i7' + g) = ±fi?g — 2t cos(A;i?a), (18) 

where the Fermi velocity is given by vp = 2tasm{kpa). The second term is a constant and 
can be shifted away by redefining the energy zero point. We will drop it throughout the rest 
of the paper. 

In one dimension, the definition of exchange statistics is ambiguous, since the no double 
occupancy condition excludes the possibility to physically exchange spatial position of two 
particles. This makes the statistics of fermionic particles lose its absolute meaning and 
make an alternative description in terms of bosons possible. This situation is quite different 
from that of the three dimensional case, where the exchange statistics of particles has an 
absolute meaning. In two dimensions, the definition of particle statistics only marginally 
makes sense and we can transmute the statistics arbitrarily by attaching the Chern-Simons 
fiux to particles (the composite of particle and fiux is dubbed as anyonSil). The statistics 
transmutation procedure in one dimension is called bosonization; it has been widely used in 
exploring the physics in one dimensional systemJ^l'B. 
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In practice, it is convenient to discuss the bosonization in real space. To do so, we expand 
the lattice fermion in terms of continuum fields 

c,. = v^[V^ij.(x)e*'=^^ + V'L.(a;)e-*'=^T (19) 



c 



i = + ^{A^)e''n. (20) 



where x = ja is used. After linearization and dropping fast varying terms, we get the low 
energy effective Hamiltonian for the Hh term as 

Hh= dxHhix) = -vfYI / dx[4j\^Jd^ipR„ - ipljd^tpLa], (21) 

describing a one-dimensional relativistic Dirac particle in continuum. 

In the following, we use the bosonization rule to bosonize the t — model: 

M^) = ^e^^-(-), (22) 
V/vra 

where t]p {P = R/L = +/—,a,...) are the Klein factors to maintain the anti-commuting 
relations between particles on different sites. We can also realize the Klein factor as Majarona 
fermions which satisfy the following anti-commuting relations 

{r]r,r]s} = '^Srs- (23) 

The introduction of l/v^27ra in Eq. (^) maintains the correct dimension for the field ipp{x) 
which has dimension [length]'^^"^ (see Eq. ([19|)). The fields are angular variables and thus 
dimensionless. To get the correct anticommutation for fermionic fields, we also require the 
bosonic fields (f)p{x) satisfy the following commutation relations 

[(f>Pa, 4>Pa'] = iPuS^^'E^x - x'); (24) 

[0ij<T, 0La'] = i-n^S^a', (25) 

where e{x) is the Heaviside step function. 

Using these bosonization rules, we get the bosonic description for the hopping term ( pi] ) 

as 
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4n 



^ J dx[{d^<f)R^y + {d^(j)Lay]. (26) 



For later convenience, let us introduce a pair of conjugate non-chiral bosonic fields for each 
species: 

(pa = (pRa + 4>La; (27) 

= - (28) 



which satisfy 



[(j)^{x),e^>{x')] = -iAn5^^,e{x - x'). (29) 



To organize the spin and charge modes more elegantly, we introduce new pairs of dual fields 
as follows: 

0c(a:) = -^(0t + (Ps{x) = ^{(p^ - (30) 

e^ix) = ^(^T + Oi); Os{x) = -^(^t - Oi). (31) 

Here the introduction of numerical factor l/\/2 is to maintain the commutation relation in 



Eq. (pQl). The subscript "s" means the spin mode and "c" the charge mode. Using these 
spin-charge separated modes, the bosonic Hamiltonian ( p^ ) can be cast into the following 
form 

Hh = ^l cia:[(9,.0e)' + (9.0,)2 + idj.y + (dj^n (32) 

To bosonize the induced interaction terms Hi 2, we note that, roughly speaking, both 
the Hi and H2 terms are of the type of Hubbard-like on-site interactions in the continuum 
limit and, therefore, provide interactions to renormalize the charge/spin velocity and the 
controlling parameters (i.e. K^/s] see below) and the cosine term in the spin sector. By 
taking microscopic details of these two terms into account, we get an extra numerical factor 
cos(A;i;'a) = cos(7rz//2). Both terms have the same bosonized form. Namely, the bosonized 
form of Hi + H2 is 
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Hi + H2 = 2Hi (33) 

where P = VriVl^VriVli, since = 1, we get V = ±1; in the following, we will take 
V = +1. 

Now we come to discuss the six-fermion term in the continuum limit; after a straight- 
forward but tedious calculation we get 

Hs= - v^tA^a^J30s(fc^a) j ^^^^^^^^ ^ _ ^^^^ 



+ 



'^'^3^^^ / dxdMx + a) cos(v^^.) 

J dxdx9s{x + a) sin(v^6's). 



tA^ sm{kFa) 



To get a sensible continuum limit, we have to take a — 0, but keep ta finite consistently. An 
elegant way to accomplish this is to use the following operator product expansion (OPE): 

dziOc{zi)dz29c{z2) ~ --, — r^; (35) 

[Zi — Z2) 
Pi 

dMz) sin[V2^,(0)] ~ cos[V29s{Q)l 

with all other OPEs being regular. We finally get the continuum limit of the six-fermion 
term to be 

H3 = I dxcos{V2e,). (36) 

The last thing in bosonizing the t — Jz model is to bosonize the magnetic interaction 
U{Jz) in Eq. ([T5|) . We decompose it into the following combinations: 

UiJz) = ^ 51 Kt - '^ii)(%+iT - ^j+n) (37) 
j 

j j 
The terms in the first bracket are Coulomb interactions between the electrons on different 

sites; in the continuum limit we get its bosonized form as 
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^ EK'T^i+iT + ^jl^j+n) = J dxlidJcY + {dJsf]. (38) 



The terms in the second bracket of eq. (|37|) are the Hubbard-hke interactions in the contin- 
uum hmit; the bosonization procedure gives 

j 

(39) 



Combining Eqs. (|38|) , (39) with Eq. (0), we find that the terms involving the charge 



variable 9c exactly cancel and we get the bosonized form of U{Jz) as 

U{J.) = ^/ dx{dj,f - -^-'"^'J,^^^^"^ / dxcos{V2e,). (40) 

It is worth noting that the absence of charge variables in the U{Jz) term is natural, since 
we are dealing with pure magnetic interactions. 

Therefore, after collecting all the results, the low energy effective Hamiltonian for the 
bosonized form of the t — model is nicely written as 

Htj^ = j dx{Hc + H,), (41) 

where the Hamiltonian for the spin sector (Hg) and the charge sector (He) are given by 

i7, = |M + i-(9A)^], (42) 

Hs = |[i^.n2 + ^^{dJs)'] + ^ cos(v^^,), (43) 

where He = ■^dx4>c and = —^dx(f)s are the conjugate momenta for the charge field 9c 
and spin field 9s respectively. The effective coupling constant is given by 



Vp 



8Acot — H cos TTZ/ — 



(44) 



2 Vp 71 

The velocities Vc/s are renormalized by magnetic interactions and the interactions induced 
by the no double occupancy conditions: 
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I J^a 4A 
Vs = vf\ 1 H 

V TVVf Tf 



Vc = Vf\ I H cot(-z/); 
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(45) 



(46) 



The controlling parameters Kc/s are given by 



An 



(47) 



V'l + ^cotdz.)' 



/! + Ail _ 4A ^.^ X 



(48) 



In passing, we would like to stress that the above results are derived for small A and J^a. 
However, the general result of a renomalization of Vc/s and K^/s, but with no other changes, 
is expected to be valid more generallySB In other words, the functional forms of the low 
energy effective Hamiltonians H^/s, being basically dictated by the symmetry requirements, 
survive even if the interactions are strong, while the above values of Vc/s and Kc/s are not 
universal. Therefore, we conclude that if we adiabatically continue the value of A to unity, 
the low energy effective Hamiltonian of the t — model should be of the same form as the 
above charge and spin Hamiltonians -ffc/s5 with the renormalized values of v^/s and K^/s not 
restricted to those given by Eqs. (^) and (^Tj). 



Now we are in the position to discuss the possible phase diagram for the t — model. 
For convenience, we only discuss the anti-ferromagnetic case, namely, we assume Jz > 0. 

At first, we notice that the spin and charge degrees of freedom are well separated just 
like what happened in other ID interacting models. However, from the expression for the 
controlling parameters Kc/s, we have already seen the interesting interplay between hopping 
and magnetic interactions. The phase diagram is determined by the competition of above 
two energy scales( t and Jz)- This is quite different from the case of the Hubbard model 
or the XX Z model where the controlling parameter is only determined by the interaction 
strength. However, the charge sector is massless, described by a quadratic Hamiltonian with 



C. The phase diagram 
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no mass term. This means that charge excitations are gapless and the charged sector of the 
system is metalhc. In contrast, the knowledge on the fate of the spin sector needs more 
work. The situation is similar to that of the Hubbard model. 

The fate of the spin sector is determined by the well-studied sine-Gordon Hamiltonian. 
In the spin sector we have the following renormalization group equations (RGE)i2lii: 

where 5 is a positive, regularization dependent parameter. With these two RG equations in 
hand, we can readily analyze the phase diagram for the spin sector. 

(i) When Kg > 47r and \ge\ < (-^ — 1)/V6, the spin sector flows to the fixed point line: 

9*e = 0, (51) 
Ks > 47r. (52) 

Thus we get the Luttinger liquid behavior for the spin sector. Following the Balents-Fisher's 
notation!, we say that the system is in the ClSl phase; here more generally a CniSn phase 
means a phase with m massless charge modes and n massless spin modes respectively. 

(ii) When parameters Kg and ge satisfy one of the following conditions: 

Ks < 47r, ge > 0; (53) 

or 

Ks > 47r, 9e<{^- 1)/ A (54) 

then the RGE flows toward go = +oo. In this case, the behavior of the system is overwhelm- 
ingly determined by the minima of the cosine term. For gg > 0, these minima are given 

by 

9s = V2in + ^)7r, (55) 
15 



but due to the angular nature of the variable ^s, we can have only two distinct ground states, 
distinguished by the even and odd values of n. This state is identified to be Peierls ordering 
of spin degrees of freedom. Due to quantum tunneling, degeneracy of the ground state is 
removed. Consequently, the excitations above either ground state are gapful. The dominant 
contributions to the mass gap come from the topological soliton excitations in the dilute 
gas approximation of solitons and anti-solitons. Therefore, in this phase, the spin sector is 
gapful, and we classify the phase of the system as a CIS'O phase. 

(iii) In contrast to the case (ii), if the parameter Kg and ge satisfy one of the following 
two conditions: 

Ks < 47r, ge < 0; (56) 

or 

Ks > An, ge < -(^ - 1)/VS, (57) 

then the RGE flows toward ge = — oo. An argument similar to that in the case (ii) gives the 
ground states determined by 

9s = V2nn. (58) 

In this state, we have a staggered expectation value for the 2— component of the spin. 
Therefore, the spin ordering is Neel-like. 

In summary, we construct the phase diagram for the t — Jz model in Figure 1. 



III. THE T - J^- Jz MODEL 



In this section, we will discuss the modified version (|^) of the t — J model. Again the 
change to make is in the magnetic interactions. In addition to the U{Jz) term discussed in 
the last section, we now add the XY part, U{J±), of the anti-ferromagenetic interactions: 

U{J^) = E(^J^j+i + S]Sl,) (59) 



J 



-L t t t t N 



2 
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Following the bosonization procedure presented in the preceding section, we get the 
bosonized form for U{J±) as 

f/(J^) = / dx cos{V2(l),) cos(V2^,) - -^[1 - cos(2A;Fa)] / dx cos{V2(ps). (60) 

Therefore, for the modified t — J model (|^), we have the bosonized low energy effective 
Hamiltonian 

n = I dxiH, + Hs), (61) 



where the Hamiltonian of the charge sector. He, is still given by Eq. (^21) , since the XY part 
of magnetic interactions only changes spin dynamics. In contrast, due to the extra U{J±), 
the Hamiltonian Hg of the spin sector has been drastically modified to 

Hs = ^-^[KsIll + ^{dJsr] (62) 



+ 



/ dxcos(V26s) - o / dxcos(V2(j)s 



+ I dxcos{V2<j)s)cosiV29s). 



Compared to the t — model, the Hamiltonian Hg now has two extra terms with the 
coupling constants and ge,/, respectively. The spin velocity (vg) and controlling parameter 
(Kg) are still the same as those in Eqs. (^5| - ^7|) . Using Eq. (160|) , the coupling constants g^f, 



and ogj, are determined to be 



= 8Jxasin2(|i/); (63) 
= 4J^a. (64) 

Due to the appearance of the interaction term gg^, which has a non-zero conformal spin, 
the dynamics for the spin sector becomes much more involved. When we use the scaling 
arguments to discuss the relevance of the interaction terms, we need to be more careful. 
We'd better use the RG flow for the Hamiltonian ( |0^ to discuss the details of the spin 
dynamics. Fortunately, up to one loop level, the RGE for a Hamiltonian like (|B^ have been 
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studied thoroughly, though in quite different contexttJlH^. The resulting RGE for the double 
cosine term oaj, is 



dl 



(65) 



Since we know 



i^s 47r 

4^ + 7f:^2' 



the double cosine term in the Hamiltonian (|62|) is always irrelevant. Of course, the action 
of RG will generate more terms, such as single cosine terms. However, the arguments of 
these single cosine terms are twice bigger and these terms are more irrelevant than the 
existing terms. Thus we can neglect them. This situation is quite different from that of the 
two coupled Luttinger liquid case^Hnl. Therefore, we only need to focus on the following 
effective Hamiltonian 

H, = ^[K,Ul + ^{dJs)'] + ^ / rfa;cos(v^^,) - ^ J cos(v^</),), (67) 

= |[n^ + [dAf] + / dxcos{(3A) - ^ I dxcosCPsk). 

where we have introduced tig = \^K^Ils and 9^ = The definitions of Ps and (3s are 



(3s = pKs, (3s = (68) 

It is now easy to observe that the low energy effective Hamiltonian possesses the following 
duality property: Namely, the Hamiltonian (|67D is invariant under the following transfor- 
mation 

(3s < — > ps, ge < — > -g<p- (69) 

Note that such a duality does not appear in the t — J^ model or in the Hubbard model. But 
it is also interesting to note that it appeared in the ID XY Z Thirring modefl and in the 
case of two coupled Luttinger hquids&H. 
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Compared with the sine-Gordon system, the symmetry of Eq. (^) is discrete, while there 
is a hidden U{1) symmetry in the sine-Gordon system, which reflects the U{1) invariance of 
its dual fermionic model (the Massive Thirring model). 

It is also easy to get the scaling dimension for the cosine terms of the field dg and its 
conjugate (pg as 

A*^^. (70) 

Therefore, one of the two cosine terms is always relevant, which is associated with the 
ordering of the 6- or (p-field. Let us discuss the following two different cases separately. 

(i) When the scaling dimension Ag < 2, the cos(v^^s) term is relevant. This case is 
similar to the t — case, and the system eventually flows to the spin-Peierls phase for 
(^e > or the Ising-Neel order for ge <^ respectively. 

(ii) When the scaling dimension A,^ < 2, the cos(v^0s) term is relevant. In this case, 
since is always negative, therefore the system flows toward the Ising-Neel phase only. 

In summary, we see that the phase diagram in Fig. 1 can only be partially accessed in the 
t — J model. The difference in the two cases reflects the fact that the duality transformation 
(|5D|) can only be reahzed in part of the parameter space, since the coupling constant g^ 
is definitely negative, while the coupling constant ge can be either negative or positive, 
depending on the interplay between t and J^. 



IV. CONCLUSIONS AND DISCUSSIONS 

In this paper, the phase diagram of the most general ID t — — Jz model is discussed 
based on bosonization and RGE. To make sense of the bosonization procedure for the inter- 
actions induced by no double occupancy constraints, we have introduced deformed Hubbard 
operators (^, which contain a deformation parameter A. While at A = 1 the no double 
occupancy constraints at each site are recovered, the case with a small positive A is acces- 
sible to perturbative RG analysis. Since the basic structure of the bosonized low energy 
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effective Hamiltonian is argued to be determined only by the symmetry requirements, the 
bosonized form of the low energy effective Hamiltonian with a small deformation parameter 
A is expected to survive the limit A ^ 1. However, we can not simply use the values of 
Vcfs and Kcjs to make precise predictions on the phase diagram, since these values are not 
reliable at A = 1. We should take the strategy in which both v^js and K^s are considered 
as phenomenological parameters. 

For the case with Jz ^ J±, the model is reduced to the so-called t — Jz model. In this 
case, the spin sector can flow to three distinct phases: the gapless phase, the spin-Peierls 
phase, and the Ising-Neel phase, depending on the range of the parameters, meanwhile the 
charge dynamics remains always gapless. In the case with > Jzc, where the Jzc represents 
the value to make ge = and Ks < 47r, the system flows to the Ising-Neel ordering in spin 
dynamics. We identify this phase as the so-called phase separation (PS) phase. For the 
case with Jz < Jzc and Ks < 4:1t, the spin sector eventually flows to the spin-Peierls phase 
which is gapful. We can identify this phase as a superconducting phase (SC). Finally, for 
Kg > Att, the spin sector flows toward a gapless phase and thus the system flows toward 
the Tomonoga-Luttinger liquid phase. Such a phase is consistent with the phase diagram 
constructed by Los Alamos group in Ref. where the authors mapped the t — Jz model 
into the ID XXZ model and construct the phase diagram from the knowledge of exact 
solutions for the Id XXZ model. This consistency also helps us to justify our proposal to 
use the deformed Hubbard operators and the continuation from the case of A ^ 1 to the 
desired case A = 1. 

In the opposite limit, namely J± ^ Jz, the modified t — J model can be reduced the the 
t — J_L model. In this case, we still have ge generally non-zero due to the no-double-occupancy 
induced interactions. Therefore, the phase diagram of the t — Jx model is expected to be 
similar to the case of the most general t — J± — J^ model. Namely, the system is generically 
gapful in the spin sector and thus can not flow toward the Tomonoga-Luttinger phase. This 
result is a little bit different from the naive speculation that the t — Jz model should be 
basically similar to the t — J model. From our study, we conclude that there are some 
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delicate differences between the two cases, since the XY part and the Ising part of the 
magnetic interactions play a different role in spin ordering. 
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FIGURES 

FIG. 1. The schematic phase diagram. The RGE flow gives the possible fate oit — Jz model 
as the spin-Pcicrs phase (S.P.), Ising-Neel pahse (I.N.), and Tomonoga-Luttinger phase (T.L.). 
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